Strongly interacting Fermi gases with density imbalance 
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We consider density-imbalanced Fermi gases of atoms in the strongly interacting, i.e. unitarity, 
regime. The Bogoliubov-de Gennes equations for a trapped superfiuid are solved. They take into 
account the finite size of the system, as well as give rise to both phase separation and Fulde- 
Ferrel-Larkin-Ovchinnikov (FFLO)-type oscillations in the order parameter. We show how radio- 
frequency spectroscopy reflects the phase separation, and can provide direct evidence of the FFLO- 
type oscillations via observing the nodes of the order parameter. 



Bardeen-Cooper-Schrieffer (BCS) pairing is behind 
various forms of superconductivity and superfluidity. A 
prerequisite for BCS pairing is the matching of the Fermi 
energies of the two pairing components (e.g. spin up and 
spin down electrons). In the case of spin-imbalanced (po- 
larized) Fermi energies, non-standard forms of pairing 
are predicted, such as Fulde-Ferrel-Larkin-Ovchinnikov 
(FFLO) pairing 111 |2j, interior gap superfluidity or 
breached pairing H H H m and phase separa- 
tion H 0|. These exotic quantum states have rel- 
evance to many fields of physics, e.g. superconduc- 
tors in a magnetic field, neutron-proton pairing in nu- 
clear matter, and color superconductivity in high den- 
sity QCD; for a recent extensive review see llj. The 
newly realized stro ngly interacting superfiuid Fermi gases 
EH El El El El El El El E3| offer a promising play- 
ground for the study of pairing and superfluidity with 
variable initial conditions — also imbalanced Fermi en- 
ergies. Indeed, first such experiments have recently been 
done, showing disappearance of superfluidity and vortices 
with increa sing spin-density imbalance |2l| . and phase 
separation |2lJ, |22j. Here we consider theoretically the 
density imbalanced Fermi gas in the unitarity regime, 
discuss the role of finite size effects, and show how phase 
separation and FFLO-type oscillations appear and can 
be observed in the RF-spectrum of the system. 

Atomic Fermi gases are confined in magnetic or op- 
tical traps and the harmonic trapping potential causes 
significant finite size effects. For instance, for the density 
imbalanced system, clear phase separation of the major- 
ity component, i.e. the component with the most atoms, 
towards the ed ges of the trap has been experimentally 
observed in [2lL l22j. Theoretically, the finite size can be 
taken into account by solving the Bogoliubov-deGennes 
(BdG) equations in the trap geometry. This has been 
done for the density imbalanced case in the BCS-limit 
(weak coupling limit) in (H 0, HI]. In the opposite 
limit where the coupling is so strong that dimers are 
formed and undergo Bose-Einstein condensation, the sys- 
tem has been described by a mean-field treatment of a 
bosonic condensate interacting with fermions in normal 
state within the local density approximation 26] . Here 
we consider the intermediate case, where the interactions 
are strong, but the pairing is still fermionic in nature. In 



the ultracold atomic Fermi gases this corresponds to the 
Feshbach-resonance, or unitarity, regime. This regime 
was considered both in 0, 111] , thereby our results pro- 
vide direct comparison to the experiments. 

We use the single-channel mean-field resonance super- 
fluidity Hamiltonian [U Eif 
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where Vt iap — \muJ^r 1 is the spherically symmetric har- 
monic trapping potential with m the atom mass and loq 
the trap oscillator frequency, is the chemical poten- 
tial for atoms in hyperfine state a, and U is the effective 
interaction strength. This describes a two-component 
((j =t, ij now two different hyperfine states of the atom) 
gas with s-wave contact interactions, where the strength 
of the interaction is tunable via a Feshbach resonance. 

Following the treatment in Ref. [H, we expand the 
field operators in eigenstates of the harmonic potential 

^( r ) = ^2i>nlm(r)Cnl m a = ^ R nl(r)Yi m (f)c n i ma 
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where Yi m (f) are the standard spherical harmonics and 
the radial part is given by 
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where r = ^Jmuj^r and L^ +1 ^ 2 is the Laguerre polyno- 
mial. This gives the Hamiltonian 
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where the single-particle energies are £ n ; ff = 
hujQ(2n + I + 3/2) — /in-, the matrix element 
F lin' = Jo°° drr 2 R n i{r)K{r)R n n{r) corresponds to 
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anomalous term that describes the Cooper pair field 
A(r) and J l nn , a = ^2 rT i^- a Jq°° dr r 2 R n i (r)n a / (r)R n n(r) is 
the Hartree interaction term. The order parameter is 
given by 

^ 2/ + 1 

A( r ) = U ^2 4?r Rni{r)R n >i{r)(c n ioic n noi), (3) 

nn'l 

and the fermion densities are 
^ 2i _|_ \ 

n<r(r) = -^-Rni(r)R n 'i(r)(cl lWa c n , l0a ). (4) 

nn'l 

These equations are solved self-consistently for fixed 
atom numbers and N± by varying the correspond- 
ing chemical potentials The Hamiltonian J2J is di- 
agonalised using the Bogoliubov transformation and the 
resulting eigenstates are used to calculate the excitation 
gap and density profiles from Eqs. © and The di- 
agonalisation gives rise to positive and negative eigenen- 
ergies E n i a , and since in general setting the particle-hole 
symmetry is broken, we need to keep all the solutions. 

We have solved the excitation gap A(r) and the density 
profiles n a (r) at zero temperature for 4570 atoms in the 
majority component (iVf) while the number of atoms in 
the minority component (iVjJ varies. We also tried finite 
temperatures but observed no significant changes to the 
picture. The parameters have been chosen for 6 Li in the 

unitarity limit with interaction parameter kpa = —16, 

h 2 k 2 

where the Fermi wave vector is given by „ F = Ep = 
Swo(6A r |) 1 / 3 . The resulting profiles are shown in Figs.Q] 
and |21 for several polarizations P = . 

The main qualitative features of the density profiles 
in Figs. ^ and [21 agree well with the experiments in 
Refs. |2lJ,|22j. Fig. ^ shows the calculated radial density. 
Experimentally, the densities are observed via imaging 
from one direction, leading to column integrated densi- 
ties. Corresponding integrated results from our calcu- 
lations are shown in Fig. [21 The density profiles show 
the phase separation into a superfluid core and a nor- 
mal fluid shell. The core corresponds to equal densities 
of the two components while at the edges the majority 
component atoms are dominating. In the column inte- 
grated picture Fig. the excess amount of the majority 
component at the edges of the trap leads to an effec- 
tive density difference also at r = 0. This is seen both 
in Fig. [3 and in the experiments 0, H2 ■ The bump 
in the density difference at the edge of the trap is an 
even more clear signature of phase separation. In [22j |. 
two regimes were observed, namely: below polarization 
P = 0.1, a coexistence regime where the density differ- 
ence does not show clear bumps (actually deviations from 
Thomas- Fermi profiles were used as the measure), and 
phase separation regime where these features are clearly 
visible. According to our results, there is no sharp tran- 
sition between these two regimes; the phase separation 
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Figure 1: The radial density profiles for the majority (f-state, 
shown in dashed line) and minority ( j-state, shown in dotted 
line) components n a (r). Solid line shows the density differ- 
ence as a function of distance from the center of the trap. The 
polarizations are P = 0.04 (upper panel), P — 0.17 (upper 
middle panel), P = 0.34 (lower middle panel), and P — 0.49 
(lower panel). 



starts immediately even for small polarizations, but the 
effect then may well be too small to observe. The ab- 
sence of a sharp transition does not exclude the possibil- 
ity of a cross-over behaviour, where the amount of phase 
separation starts to grow faster after a certain thresh- 
old. To this extent, we plot in Fig. 01 the excitation gap 
A at center of the trap, as well as the size (volume) 
of the superfluid core, as functions of the polarization. 
There might be a change of slopes around the polariza- 
tion P — 0.2, but these results alone are not sufficient for 
any conclusive statements about a crossover behaviour. 
However, they tell clearly that the value of the gap at 
the center tends to stay quite constant, and the effect 
of the polarization is mainly to decrease the superfluid 
core size. Indeed, the superfluid core size becomes neg- 
ligible at some polarization P between 0.6-0.9, agreeing 
well with the experiments [2]]], where P = 0.7 was found 
to be a threshold for the disappearance of the condensate 
at the unitarity limit. 

In Fig. |S1 we show how the phase separation is reflected 
in the RF-spectrum of the gas. RF-spectroscopy |3(| 
has been used ej? to observe the excitation gap of the 
system We calculate the spectra using the 

method presented in [33j. The results of Fig. [S] show a 
broad peak at finite RF-field detuning, corresponding to 




Figure 2: The column integrated density profiles for the ma- 
jority (dashed) and minority (dotted) components. Solid line 
shows the density difference as a function of distance from 
the center of the trap. The polarizations are the same as in 

Fig. El 
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Figure 3: The radial density (dashed for f-state and dotted 
for j-state) and the gap (solid line) profiles for polarizations 
P = 0.34 (left) and P = 0.88 (right). 

paired atoms, and the detuning is related to the pairing 
energy A. Note that there are equal amounts of paired 
atoms for both components. The narrow peak at zero de- 
tuning corresponds to the nonpaired majority component 
atoms at the edge of the trap. This could be a probe of 
phase separation, complementary to the straightforward 



Figure 4: The value of the order parameter at the center 
of the trap (squares) and the fraction of the superfluid core 
as compared to the Fermi sphere of the major component 
(circles) plotted as functions of polarization P. 




Figure 5: The RF-spectra for different components. Left 
panel shows the spectra for both majority (dotted) and mi- 
nority (solid) components for polarization P — 0.34 and the 
right panel shows the spectrum for the minority component 
for polarization P = 0.88. The calculations are done at T = 



observation of the density profiles, since it does not suf- 
fer from the effects of column integration and provides 
a direct comparison between the amounts of the paired 
and nonpaired atoms. 

The gap and density profiles in Fig. [3| show the os- 
cillation of the order parameter and the density as a 
function of the radial coordinate. Such oscillations have 
been accounted 00113 for a FFLO-type phase. The 
FFLO-state in homogenous space leads to oscillations of 
the order parameter, and is by definition pairing with 
unequal number of particles. In homogenous space, the 
FFLO pairing starts only after a critical polarization. In 
our results, oscillations are also visible for small polar- 
izations, although as tiny effects. This is understandable 
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in the sense that, as the trap favours phase separation, 
the local polarization at the edges of the trap becomes 
very easily of considerable size. Therefore locally one can 
fulfill the FFLO condition of exceeding a critical polar- 
ization, even when the total polarization is small. One 
could interpret the results in the following way: the trap 
tends to enforce a normal BCS state at the center of the 
trap and a FFLO-type state at the edges, and the signif- 
icance of the latter grows with the total polarization. Is 
this FFLO-type state observable? It may have existed in 
the experiments [2llE^. The oscillations of the order pa- 
rameter are accompanied by oscillations of the densities, 
and therefore, in principle, one could observe the FFLO 
characteristics from the density profiles. The column in- 
tegration, however, tends to wash out the oscillations as 
can be seen by comparing Fig.^and Fig. [21 Thus, experi- 
mentally it may be difficult (although column integration 
can in principle be avoided by more advanced techniques) 
to detect the FFLO-phase from the density profiles. The 
nodes of the order parameter, however, should be visi- 
ble in the RF-spectrum of the minority component; they 
produce a peak at zero detuning, reflecting a finite num- 
ber of non-paired atoms also in the minority component. 
This too, may be a small effect for some parameters, 
e.g. in the left panel of Fig. such zero detuning peak 
is not visible in the minority component. However, for 
parameter values that also produce a more prominent 
oscillation of the order parameter (Fig. [3] right panel), 
the zero detuning peak becomes clearly visible, see the 
right panel in Fig. 03 This is a direct evidence for the 
nodes of the order parameter. Situations where such sig- 
natures are large enough to be observed can probably be 
achieved experimentally. For instance, we were restricted 
to spherical geometry due to computational reasons, but 
a cigar-shaped system is likely to display more prominent 
oscillations. 

In summary, we considered trapped, strongly interact- 
ing Fermi gases with unequal populations of the pair- 
ing components. We relate our findings to recent ex- 
periments and suggest new ways of observing the phase 
separation and, especially, FFLO features. The system 
seems to be suited for detailed studies of exotic forms 
of fermion pairing. Our results show that the trapping 
potential affects the system in an essential way; spatial 
regions with different pairing characteristics tend to form 
and finite size effects have to be carefully taken into ac- 
count in understanding the system. 
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